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One-half  crack  length 
Young's  Modulus 

Stress  intensity  factor  =  aAc" 

Strain  hardening  exponent 

Cylindrical  coordinates  measured  from  crack  tip 

Secant  modulus  =  ae/ey 

Displacement  vector 

Strain  energy  density 

Distance  in  front  of  deformed  crack  tip 

Distance  perpendicular  to  deformed  crack  tip 

Strain  tensor 

Plastic  strain  tensor 

Initial  yield  strain 

Effective  plastic  strain 

Effective  total  strain 

Deformed  crack  root  radius 

Poisson's  ratio 

Stress  tensor 

Hoop  stress 

Effective  stress 

Tensile  yield  stress 

Remote  tensile  stress 

Airy  stress  function 
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INTRODUCTION 


In  recent  years  Cherepanov  [1],  Rice,  [2,3],  Hutchinson  [4,5],  and 
Rice  and  Rosengren  [6]  have  shown  the  asymptotic  behavior  of  stress  and 
strain  fields  surrounding  sharp  crack  tips  in  plane  strain.  Using  these 
studies  as  a  guide,  full  field  solutions  with  finite  elements  have  been 
obtained  by  Levy,  Marcal,  Ostergren  and  Rice  [7]  and  Hilton  and  Hutchinson 
[8].  These  two  studies  give  accurate  near  and  far  field  behavior  due  to 
the  inclusion  of  singular  elements  reflecting  plasticity  at  the  crack 
tip.  Other  numerical  solutions  by  Marcal  and  King  [9],  Mendel  son  [10], 
Swedlow  and  coworkers  [1 1 ,1 2]  and  Tuba  [13]  show  qualitative  features  of 
the  near  field,  but  may  not  yield  accurate  stress  field  definition  due 
to  the  large  gradients  there. 

In  order  to  have  an  accurate  description  of  the  near  field  surrounding 
crack  tips,  and  hence  a  good  understanding  of  the  mechanisms  of  failure. 
Rice  and  Johnson  [14]  have  pointed  out  that  crack  tip  blunting  must  also 
be  included.  Their  analysis  accounted  in  an  approximate  manner  for  the 
blunting  at  the  crack  tip  and  for  strain  hardening  in  the  plastic  zone. 

As  a  result  they  showed  that  the  stresses  near  the  crack  tip  are  indeed 
finite  and  that  the  maximum  hoop  stress  occurred  at  some  small  distance 
from  the  deformed  crack  tip.  A  finite  deformation  analysis  by  McGowan 
and  Smith  [15]  of  blunted  cracks  in  a  linear  (stress-strain)  incompressible 
material  shows  the  same  general  behavior.  The  maximum  hoop  stress  occurs 
in  front  of  the  blunted  crack  tip  and  the  magnitude  is  independent  of  the 
remote  loading. 

The  purpose  of  the  present  study  is  to  gain  a  full  field  solution 
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around  a  blunted  crack  tip  in  a  strain  hardening  incompressible  material 
under  Mode  I  loading.  This  work  will  provide  an  accurate  description  of 
the  stress  and  deformation  fields  immediately  surrounding  the  blunted 
tip,  and  thereby  gain  insight  to  fracture  behavior.  Deformation  theory 
of  plasticity  with  a  Mises  yield  condition  is  used.  The  resulting  set 
of  equations  is  solved  for  the  blunted  crack  tip  in  the  deformed  state 
under  load  by  finite  differences.  The  linear  theory  of  Inglis  [16]  gives 
the  necessary  asymptotic  boundary  conditions. 

An  initial  goal  of  the  present  study  was  to  gain  a  more  complete 
understanding  of  the  near  field  behavior  of  stress  freezing  photoelastic 
materials  above  critical  temperature;  however,  this  study  should  give 
considerable  insight  to  the  general  behavior  of  engineering  materials 
under  Mode  I  loading. 

FORMULATION  OF  THE  PROBLEM 

For  small  strains  the  compatibility  equation  in  two  dimensions  is 

ell,22  +  e22,ll  “  2el 2,12  =  0  ^ 

Using  a  deformation  theory  of  plasticity,  the  constitutive  equations 
are  [10] 

ell  =  F  (all  "  y(°22  +  a33^  +  ^  (all  "  1/2(a22  +  a33^ 

e22  =  F  (a22  "  y(all  +  a33)}  +  ^  (a22  '  1/2(all  +  a33})  (2) 

e33  =  f  (a33  -  ^(°n  +  a22))  +  (a33  “  1/2(all  +  a22^ 
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where  eT  = 

T  p  E 

Introducing  equations  (4)  into  the  compatibility  equations  (1)  they 
become 


T  (a11,22  +  °22,11  "  all,ll  "  a22 ,22  "  4a12,12* 

+  2  (y)>2  ^ al 1 ,2  “  °22,2  "  2a12,l^  +  2  V,1  ^22,1  ‘  CT11,1  ^ 

"  2a1 2,2^  +  ^11  “  °22^M,22  "  M,ll^  "  4°1 2  ^  ,12  ~  ° 

where  T  =  secant  modulus  =  ae/e-p. 

Introducing  the  Airy  stress  function  to  satisfy  equilibrium 

CT11  =  $,22  ;  °22  =  $,11  ;  al 2  =  _$,12  (6) 

Substituting  equation  (6)  into  equation  (5),  the  governing  equation  for 
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the  field  becomes 


f  ^ ,2222  +  *,1111  +  24>,1 1 22)  +  2(T)  ,1  (4>,m  +  *,122* 

+  2  ^,2  (*,222  +  *,112*  +  ^*,22  "  $,11)[(T),22  '  ,11] 


+  4  ^T* ,12  *,12  0 


with 


VTTt(),ll  -‘,22>2  +  4*,,223V2 

For  this  study  a  Ramberg-Osgood  material  will  be  used 


(7) 


(8) 


§“  eT  =  ae/oo  +  bat(Vao)1/n  '  1] 

°  (9) 

or  j  =  1  +  ba[(ae/a0)^1_n)/n  -  aQ/ae] 

where  b  =  0  if  ag  <  aQ 
b  =  1  if  ae  >  a0 

Thus,  the  governing  equation  (7)  will  be  solved  subject  to  the  consti¬ 
tutive  laws  (equations  (8)  and  (9)). 

The  geometry  of  the  blunted  cracks  in  the  deformed  state  under  Mode 
I  loading  will  resemble  small  elliptical  perforations  as  shown  in  Figure 
1.  The  size  of  the  deformed  crack  tip  root  radius  will  be  determined 
through  integration  of  the  strain  displacement  relationships 


w.  .  +  w.  .  =  2e. • 
1 3 


(10) 


The  affected  strain  hardening  region  will  be  divided  into  a  small 
grid  utilizing  elliptical  coordinates  and  the  governing  set  of  equations 
will  be  solved  through  the  method  of  finite  differences.  At  some  dis¬ 
tance  from  the  deformed  crack  tip  the  linear  solution  of  Inglis  [16]  will 
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apply.  The  stress  at  the  outer  boundary  of  the  inner  strain  hardening 
region  will  be  then  matched  to  the  Inglis  solution.  The  outer  boundary 
will  be  enlarged  until  there  is  no  change  in  the  inner  stress  field.  (A 
detailed  description  of  the  solution  procedure  is  given  in  [17].) 

PRESENTATION  OF  RESULTS 

The  stress  and  displacement  fields  in  the  field  surrounding  a  de¬ 
formed  crack  tip  in  a  strain  hardening  material  which  is  incompressible 
in  both  the  elastic  and  plastic  regions  are  examined.  Strain  hardening 
exponents  of  0.2  through  0.01  are  presented.  The  range  of  initial  yield 
strain  values  is  from  0.01  through  0.0001.  The  value  of  a  in  the  effec¬ 
tive  stress-effective  strain  relationship,  equation  (9),  is  taken  to  be 
1.0  in  this  study.  (The  authors  have  found  that  small  changes  in  a  and 
y  do  not  influence  the  solution  significantly.)  The  "linear"  results 
reported  here  are  those  of  Inglis  [16]  for  a  deformed  crack  tip  in  a 
linear  material.  The  "singular"  results  are  those  corresponding  to  a 
crack  which  has  no  root  radius  in  a  linear  material. 

The  plastic  zone  shape  for  the  smallest  ellipse  investigated  (p  = 
0.0018  (K /aQ)2  )  is  shown  in  Figure  2.  Note  that  with  decreasing  hardening 
(n  0)  the  plastic  zone  grows  in  maximum  extent  and  leans  progressively 
in  the  direction  of  crack  propagation.  For  comparison,  the  singular 
plastic  zone  from  McClintock  and  Irwin  [18]  and  that  of  Levy  et  al  [7] 
for  a  non-hardening  (n  =  0)  material  are  also  shown.  As  shown  by  Figure  2, 
the  plastic  zone  shape  predicted  by  McClintock  and  Irwin  [18]  is  approached 
by  the  present  study  as  n  <*>.  The  difference  between  the  plastic  zone 
shape  predicted  by  Levy  et  al  [7]  and  the  present  study  for  n  =  0.01  is 
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primarily  due  to  the  inclusion  of  blunting  effects  in  the  latter;  the 

difference  should  be  negligible  as  eQ  -*•  0. 

The  effective  stress  (ag)  is  shown  versus  the  distance  ahead  of  the 

deformed  crack  tip  in  Figure  3  and  the  effective  stress  is  shown  versus 

distance  perpendicular  to  the  deformed  crack  tip  in  Figure  4.  These 

figures  show  that  in  both  cases  the  effective  stress  varies  as  ■■  1  +y  in 

r  ' 

the  plastic  zone.  It  can  be  shown  that  the  strain  energy  has  the  form: 


Therefore,  the  strain  energy  varies  approximately  as  1/r  in  the  plastic 
zone.  This  was  a  key  assumption  in  the  analysis  of  Rice  and  Rosengren  [6] 
and  Hutchinson  [4]. 

The  hoop  stress  (cryy)  in  front  of  the  crack  tip  is  shown  for  various 
values  of  yield  strain  for  n  =  0.2  in  Figure  5,  for  n  =  0.1  in  Figure  6, 
and  for  n  =  0.01  in  Figure  7.  As  shown  in  these  three  figures,  the 
hoop  stress  is  substantially  reduced  near  the  crack  tip  because  of  blunting 
and  strain  hardening;  with  the  maximum  hoop  stress  developed  forward  of 
the  crack  tip.  The  analysis  of  Rice  and  Johnson  [14]  gives  the  same 
qualitative  behavior;  the  correlation  is  believed  to  be  quite  reasonable 
in  view  of  the  several  approximations  involved.  For  a  non-hardening 
material  Rice  [2]  has  shown  that  the  maximum  hoop  stress  is  2.97  o0.  The 
hoop  stress,  as  predicted  by  Levy  et  al  [7],  approaches  this  limit  at  the 
crack  tip  as  shown  in  Figure  7.  The  hoop  stress  distribution  of  the 
present  study  in  Figure  7  reflects  the  presence  of  blunting  and  should 
coincide  with  the  work  of  Levy  et  al  [7]  as  eQ  ->  0. 

Figure  8  shows  the  variation  of  maximum  hoop  stress  with  initial 
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yield  strain  for  varying  hardening.  As  shown  in  the  figure,  blunting  alone 
(linear)  makes  the  peak  hoop  stress  finite  and  the  inclusion  of  finite  defor¬ 
mations  [15]  reduces  the  magnitude  somewhat.  However,  the  effects  of  blunting 
and  plasticity  taken  together  are  significant:  the  peak  hoop  stress 
is  reduced  by  a  factor  of  10  from  that  with  blunting  alone.  From  Figure 
8,  one  observes  the  peak  hoop  stress  to  be  3aQ  to  7 aQ  depending  upon  n 
and  oq/E.  The  peak  hoop  stress  increases  with  r  and  decreases  with  aj E. 

(The  large  value  of  peak  hoop  stress  compared  to  the  uniaxial  yield 

stress,  a  ,  is  believed  due  to  the  presence  of  triaxiality  in  the  crack 
o 

tip  region.) 

The  crack  tip  displacement  in  the  direction  of  propagation  (which  is 
also  the  deformed  crack  root  radius,  p)  is  shown  in  Figure  9  for  varying 
initial  yield  strain  and  hardening  exponent.  The  present  study  predicts 
that  the  forward  crack  tip  displacement  increases  with  a0/E  and  decreases 
with  n.  For  comparison  one-half  the  crack  tip  opening  displacement  cal¬ 
culated  by  Levy  et  al  [7]  is  shown.  The  forward  crack  tip  displacement 
as  predicted  by  the  present  study  and  the  work  of  Levy  et  al  [7]  show 
parallel  behavior,  although  they  are  separated  by  some  distance.  This 
disagreement  is  believed  due  to  the  shape  of  the  crack  tip  being  ellipti¬ 
cal  in  the  present  study  instead  of  cylindrical. 

Included  also  on  Figure  9  is  the  width  of  the  "transition"  or 
"stretch"  zone  which  exists  on  the  fracture  surface  between  the  cracked 
and  the  overload  regions  in  fatigue.  As  Broek  [19]  has  discussed,  the 
depth  of  this  transition  zone  is  the  crack  tip  opening  displacement, 
and,  therefore  the  width  is  the  tip  forward  displacement. 

Examination  of  the  figure  shows  the  correlation  between  the  forward 
tip  displacement  and  failure.  The  measurements  of  the  stretch  zone 
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falls  close  to  n  =  0.2.  For  the  steels  and  aluminums  shown  values  of  n 
around  0.05  have  been  reported  in  [19],  [20],  and  [21].  However,  it  is 
known  that  for  this  class  of  materials  the  value  of  n  varies  with  plastic 
strain  [22].  For  large  plastic  strain  (e  >  10%),  the  strain  hardening 
exponent  is  close  to  0.2  as  shown  by  Jones  and  Brown  [23]  for  4340  steel. 
The  strains  in  the  tip  region  are  clearly  greater  than  10%  so  that  the 
agreement  between  the  measurements  and  the  analysis  appears  quite  rea¬ 
sonable.  The  scatter  band  shown  on  the  figure  is  an  indication  of  the 
span  of  actual  measurements  (authors  typically  report  a  40%  variation). 


DISCUSSION 

Study  of  the  stress  fields  surrounding  deformed  crack  tips  indicate 
that  the  stresses  in  the  plastic  zone  will  approach  the  l/rn//n+1  varia¬ 
tion  as  eQ  -*■  0  as  determined  in  the  analysis  of  Hutchinson  [4]  and  Rice 
and  Rosengren  [6].  However,  as  mentioned  in  the  work  of  Levy  et  al  [7], 
the  J-integral  may  not  be  path  independent  in  the  plastic  zone;  and 
therefore,  a  full  field  solution  is  desirable  in  order  to  evaluate  the 
strength  of  the  stress  variation  in  the  plastic  zone.  For  the  non-hardening 
case,  Levy  et  al  [7]  and  subsequently  Rice  and  Tracey  [24]  did  a  full 
field  solution  of  a  sharp  crack  tip  utilizing  singular  elements.  The 
differences  in  the  present  study  and  [7]  and  [24]  is  primarily  due  to 
blunting,  and  partially  due  to  use  here  of  0.5  as  an  elastic  Poisson's 
ratio.  These  differences  should  vanish  as  eQ  ->-  0. 

The  analysis  of  Rice  and  Johnson  [14]  was  an  attempt  to  account  for 
the  effect  of  blunting  near  the  crack  tip  and  used  the  results  of  [4] 
and  [6]  as  boundary  conditions.  Although  approximate,  the  analysis  of 
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Rice  and  Johnson  [14]  is  in  qualitative  agreement  with  the  present  study 
and  the  stresses  near  the  crack  tip  are  substantially  reduced.  Also  their 
analysis  does  predict  the  maximum  hoop  stress  occurring  in  front  of  the 
crack  tip. 

Previously  McGowan  and  Smith  [15]  performed  a  finite  deformation 
analysis  of  the  region  surrounding  deformed  crack  tips  for  a  linear 
(stress-strain)  material.  The  results  of  the  finite  deformation  work 
showed  that  the  maximum  hoop  stress  occurred  in  front  of  the  deformed 
crack  tip  also.  It  was  determined  that  the  stress  distribution  around 
the  crack  tip  was  "similar";  in  the  sense  that  one  stress  distribution 
could  be  used  to  describe  the  response  of  the  material  under  load.  The 
size  of  the  affected  zone  would  depend  upon  the  load  and  crack  length  through 
K.  The  self-similarity  of  the  stress  field  was  a  direct  result  of 
the  blunting  process,  and  would  be  expected  to  remain  as  long  as  the 
affected  zone  stayed  small  with  respect  to  the  crack  length,  thickness, 
or  any  other  in-plane  dimension. 

The  behavior  is  quite  similar  for  a  power  hardening  material.  The 
stress  field  is  self-similar  with  the  size  of  the  affected  zone  varying 
with  K.  The  maximum  hoop  stress  will  only  then  be  a  function  of  the  ma¬ 
terial  properties  E,  n,  and  a^.  The  stretching  of  the  similar  stress 
distribution  will  depend  upon  K  as  well  as  the  other  ma¬ 
terial  properties.  One  may  conjecture  that  failure  would  depend  upon  the 
growth  in  size  of  a  critical  dimension  such  as  plastic  zone  size  which 
increases  with  K. 

Wells  [25]  and  others  have  used  the  crack  opening  displacement 
as  a  fracture  criterion.  Broek  [19]  has  used  this  concept  to  correlate 
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the  depth  of  transition  zones  in  aluminum  with  fracture  toughness. 

The  present  study  shows  good  correlation  of  fracture  toughness  and 
transition  zone  width.  Krafft  [26],  Hahn  and  Rosenfield  [27]  and  Rice 
and  Johnson  [14]  have  all  shown  good  correlation  of  plane  strain  frac¬ 
ture  toughness  with  some  minute  particle  size  or  process  zone  size  for 
specific  cases. 

SUMMARY  AND  RECOMMENDATIONS 

Following  the  pioneering  studies  of  Hutchinson  [4],  Rice  and  Rosen- 
gren  [6],  Levy  et  al  [7]  and  Hilton  and  Hutchinson  [8],  the  authors  have 
obtained  a  full  field  deformation  plasticity  finite  difference  solution 
to  the  Mode  I  plane  strain  problem  including  the  effects  of  blunting. 

The  material  was  incompressible  in  both  the  elastic  and  plastic  regions, 
and  followed  a  power  hardening  rule.  Stress  and  displacement  fields 
surrounding  the  deformed  crack  tip  are  presented,  and  are  found  to  com¬ 
pare  favorably  both  with  the  analysis  of  other  investigators  as  well  as 
experimental  results.  Because  of  the  improved  accuracy  expected  from  a 
full  field  solution,  it  would  be  appropriate  to  incorporate  such  a  solu¬ 
tion  into  theories  concerning  void  coalescence  and  final  instability 
(such  as  the  landmark  studies  of  McClintock  [28],  Rice  and  Johnson  [14] 
and  Rice  and  Tracey  [29]).  Efforts  are  currently  being  devoted  to  such 
an  approach. 
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Figure  2  Plastic  zone  shape  for  the  blunted  crack  for 
several  strain  hardening  exponents  with  a 
deformed  root  radius  of  p=Q.Q018(K/oo)2. 
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Figure  3  Effective  stress  distribution  forward  of  the 
blunted  crack  tip  for  several  strain  hardening 
exponents  with  a  crack  root  radius  of  p=0.0018 


Figure  4  Effective  stress  distribution  perpendicular  to 

the  blunted  crack  tip  for  several  strain  hardening 
exponents  with  a  crack  root  radius  of  p=0.0018 
(K/%P. 
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Figure  7  Hoop  stress  distribution  forward  of  the  blunted  crack  tip  for  several  values  of  initial 
yield  strain  with  a  strain  hardening  exponent  of  n=0.01. 


Variation  of  hoop  stress  maximum  with  initial  yield  strain  for 
several  values  of  strain  hardening  exponent. 
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Figure  9  Forward  crack  tip  displacement  variation  with  initial  yield 
strain  for  several  values  of  strain  hardening  exponent. 

Also  shown  are  experimental  measurements  of  'stretch'  zones 
in  various  steels  and  aluminums. 


